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Abstract 

We study gauge theory operators which take the form of a product of a trace 
with a Schur polynomial, and their string theory duals. These states represent 
strings excited on bubbling AdS geometries which are dual to the Schur polyno- 
mials. These geometries generically take the form of multiple annuli in the phase 
space plane. We study the coherent state wavefunction of the lattice, which labels 
the trace part of the operator, for a general Young tableau and their dual de- 
scription on the droplet plane with a general concentric ring pattern. In addition 
we identify a density matrix over the coherent states on all the geometries within 
a fixed constraint. This density matrix may be used to calculate the entropy of a 
given ensemble of operators. We finally recover the BMN string spectrum along 
the geodesic near any circle from the ansatz of the coherent state wavefunction. 
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1 Introduction 



The AdS/CFT correspondence is an invaluable tool [T], not only enabling us to study 
the strong coupling limits of a wide variety of gauge theories, but also allowing us to 
study deep questions arising in the study of quantum gravity. 

Certain heavy supersymmetric states in the gauge theory can be dual to geometric 
backgrounds in string theory. An example is the set of 1/2 BPS chiral primary states 
built with one complex scalar [2] , [3] , [1] , [9] , [in] . These operators can be labelled by 
Young tableaux or Schur polynomials which themselves form an orthonormal basis. 
When the charges of these operators are very large, they become dual to a variety of 
geometries, which have bubbling type boundary conditions in some loci of the spacetime 
far from the boundary, that encode the mapping from the gauge theory states. Besides 
the 1/2 BPS chiral primary operators, there are also other types of operators with 
corresponding bubbling geometries, e.g. [H] , [12] , [13] , [11] in A/" = 4 SYM. 

A natural question which arises is what is the description of the non-BPS excitations 
on top of the above states in both the gauge and string theory sides? On the gauge 
side, we study operators which are products between a Schur polynomial and a single 
trace operator. We may view the Schur polynomial as the background geometry and 
the single trace part as a further stringy excitation on the geometry. On the string side, 
we study closed strings excited on the dual geometries, described by a general Schur 
polynomial or Young tableau. In particular we study the backgrounds with general 
concentric ring distributions [1] in the phase space plane. The size of the S'^'s in the 
LLM metric are determined by the position in the droplet distribution, in which the 
edges correspond to the vanishing of one or the other 5"^. 

One way to study such non-BPS states is to look at the dilatation generator acting 
on them. The dilatation operator acting on the non-BPS state is naturally described 
in terms of the hamiltonian on a lattice, which itself can be described in a variety of 
limits. The bosonic occupation numbers of the lattice sites correspond to the number 
of complex scalar fields in the trace part of the operator in between the impurity 
fields. This has been studied before in the context of certain background geometries 
[5] , [6] , [7] , [8] . We will be particularly interested in studying the coherent states of such 
a lattice, though for more general background geometries. In the continuum limit of 
the lattice we obtain a sigma model that matches with the sigma model of the closed 
string living in the dual geometry. We can also study the lattice with a small number 
of sites, corresponding to short BMN type strings. 

In [6] an analysis was made of a reduced subspace of possible diagrams. In this work 
the authors analyzed the coherent state on a disk plus ring distribution and proved 
that the normalized wavefunction constructed from this coherent state knew about the 
geometry and topology of the string theory background. In the current work we fully 
generalize this description to include any geometries which can be described by axially 
symmetric solutions in the phase space plane. 

Furthermore we will look at the entropy of particular ensembles of such states. In 
the 1/2-BPS limit all non-equivalent geometries are orthogonal to one another, meaning 
that an ensemble made of such solutions does not mix. The breaking of supersymmetry 
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via the stringy excitation/trace operator allows for a mixing between elements of the 
gravitational ensemble and a non-diagonal reduced density matrix. 

The organization of this paper is as follows. In section 2.1, we introduce the non- 
BPS operators that we will study, the lattice labelling, and the associated oscillator 
algebra. In section 2.2, we analyze the coherent state wavefunction of the system, the 
continuum limit of the lattice and sigma models, and also the properties of the coherent 
state wavefunction, such as the norm and the average occupation number. In section 
2.3, we describe a density matrix over these coherent states with a constrained set of 
Young tableaux, and from this we are able to calculate the entropy of such an ensemble 
of states. In section 2.4, we derive the BMN states from the ansatz of the coherent 
states and show that in the appropriate limits we recover the results of [Ej. Finally, in 
section 3, we briefly discuss our results and conclude. We also include four appendices 
providing concrete computational details on the calculations discussed in the body of 
the text. 

2 Strings on Bubbling Geometries 
2.1 Operators 

The 1/2 BPS chiral primary operators of A/" = 4 super Yang Mills can be labelled by 
Schur Polynomials Xr{Z)i where R is the representation in which the operator is being 
traced, and Z is one of the complex scalar fields. Equivalently they can be represented 
by Young tableaux which are dual to 1/2 BPS geometries. These can be described using 
the LLM prescription as a series of concentric black and white annuli distributions. The 
Young tableau corresponding to a particular chiral primary operator may be labelled 
by the positions of its corners. We label outward corners with a "6" and inward corners 
with an "a" , going from bottom left to top right. A given diagram will come with labels 
bi,ai,b2,a2, ...,bm,cim where m is the total number of inward, or equivalently outward 
pointing corners. Such operators can be mapped into a distribution of m black rings 
where the inner and outer edges of a ring/ annulus correspond on an outward pointing 
or inward pointing edge of a Young tableau, respectively. We label these inner and 
outer radii as R2k-i and i?2fc where k labels the ring (we will use k and / to label ring 
numbers throughout). Due to quantization of the phase space area, the radius squared 
i?2fc_i of each inner circle of the black rings is proportional to the distance Cb^ from the 
lower-left point of the Young tableau to the outward corner bk, of the Young tableau, 
counting along the edge. Similarly, the radius squared i?2fc of ^^^h outer circle of the 
black rings is proportional to the distance Ca,, from the lower-left point of the Young 
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tableau to the inward corner B 



- ^2.-^, (2.1) 



k 

a, = C^,-N,, a, = J](M,+iV,). (2.2) 

i=i 

This provides a simple mapping between a Young tableau labelling and a phase space 
distribution. 




Figure 1: An example of the Young tableau description of a Schur polynomial operator in 
the gauge theory (with small N and M as an example). The numbers of rows and columns on 
each step (M^ and N^) are labelled along with the positions of the inner and outer edges of 
the associated annulus distribution in the LLM plane. The rings are labelled /c = 1 to 3 from 
the inner ring to the outer ring. The dark regions in the LLM plane are those with boundary 
value — ^ but are shaded in grey for clarity. The ring radii squared are: (1, 5, 9, 13, 15, 16). 

It should be noted that such Young tableaux do not exhaust all possible concentric 
ring distributions. Such a description will not include phase space distributions which 
start with a black diskl. These are labelled separately with a thin vertical line running 
from the bottom left hand corner of the Young tableau, corresponding to having no 
white region in the centre. The labelling of such states follows logically. The bulk 

*The coordinates of the z-plane in this paper are rescaled from the coordinates of the z-plane in 
g] by a factor, such that ^^^jW^ = ^^-^§r^- 

^ These would correspond to Mi = 0, i?i =0 limit, which will be discussed more in Appendix \K\ 
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of this paper is motivated by studying closed string states in backgrounds given by 
general axially symmetric ring distributions. 

Having described this subspace of the half-BPS sector we would like to look at 
products of a trace multiplied by an operator described by the Schur polynomials 
discussed above. Such products are in general non-BPS and we will utilize the breaking 
of supersymmetry in what follows. We label such products as: 

0{n,} = tr y^c.,Z:^ Xh{Z) , (2.3) 

where Xr{Z) is a Schur polynomial corresponding to a Young tableau of representation 
R. Wq, is a letter corresponding to impurities coming from fields which are not in the 
1/2 BPS chiral primary sector. Za is a field creating or removing a box with weight 
Caj. or Cb^. at the Young tableau corner, where the subscript aj denotes the corner ak 
or bk (see also [7]). Such an operator corresponds to a closed string, represented by the 
trace, in a background represented by the Schur polynomial. The presence of the string 
in the background leads to SUSY breaking. The operators that are eigenstates of the 
dilatation operator are superpositions in the space spanned by the basis (12. 3p . (In [6], 
product of determinants were used in place of the Schur polynomial part.) The total 
number of Zs in the trace part of (12.31) is not fixed when we mix these operators to 
obtain the eigenstates of the dilatation generator, so it is convenient to view the trace 
part of (12.31) as a lattice of L sites each with an occupation number rij. For simplicity, 
we first take Wq, to come purely from the complex scalar Y. Each operator can thus 
be represented by a state \ni,n2, ■■■.til) . The positive occupation number corresponds 
to positive angular momentum of the string along the outer circle of a black ring, while 
the negative occupation number corresponds to the negative momentum of the string 
along the inner circle of the black ring. Such a basis can be used to describe a wide 
variety of string states in 1/2-BPS backgrounds. 

The effective Hamiltonian of the lattice is derived from the dilatation generator 
A = Ao + AAi + -^"^n, where 

n>2 

Ao + AAi = tr{Zdz + Ydy) + 4^tr[Z, Y] [dz, dy] , (2.4) 



for the operators under consideration, and here A = We have 

L 

^0 = {J,, + ^nj) + L = J + L, (2.5) 
\D, = ^tr[Z,Y][dz,dy] = A-J-L = H, (2.6) 

OTT 

where in the last line we only consider the one-loop dilation generator. The Hamiltonian 
when acting on the lattice takes the form 

L 

H = E-J-L = \ ^(aj - aj+i)(a,- - a,+i) , (2.7) 
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where aj are shift operators which reduce the bosonic occupation number by one. We 
use J" to denote the total ?7(l)-charge associated with Z, and use J to denote that 
t/(l)-charge from the trace part, i.e. 

L 

J = Jtrace = ^j- (2-8) 

For a general Young Tableau, we have the following algebra 



a\-n,bk) = -^^7=^ |-n - l,6fc) , n ^ 0, (2.9) 

a\l,ai) = y yi ^ |0)fc , a\n,ak) = —r^\n-l,ak) , n > 1 (2.10) 

V N 




a^\n,ak) = ^^\n + l,ak), n ^ 0, |0, 6^) = |0, a^) = |0),, , (2.11) 

a^\-n,bk) = ^^^\-n+l,bk), n>0. (2.12) 

The state \—n, b^) corresponds to the —n^^ excitation on the inner edge of the k^^ ring 
while |n, Ofe) corresponds to the n^^ excitation on the outer edge of the k^^ ring. As 
noted before, the Cb,. and Ca,. are the distances (number of boxes) between the k^^ 
outward pointing or inward pointing corner and the bottom-left point of the Young 
tableau. The Cb^ and Ca^. are also the weights of the boxes of these Young tableau 
corners, and are the coefficients that appear in the two-point correlators for the fields 
2a j [Z]- The action of the annihilation/creation operator on different types of Za- field 
excitations differ by the coefficients ^/CbjN ■, or a/ Ca^ / N, and is otherwise similar. 
The action of aj (a^) removes (adds) a box in the appropriate corner of the Young 
tableau. More discussions of their properties are included in Appendix |X1 

In principle one can take several limits in order to study superpositions of operators 
of the form (12. 3p . and their string theory duals. The limits of interest (both on the 
gauge and string side) include: 

1. The continuous limit /sigma model limit: 

L > 1, J unfixed; A < 1, v.s. \ <^ I (2.13) 

2. The BMN/plane-wave type limit: 

L = 2, |J|>1, ^T^l; v.s. i<l (2.14) 

A 



3. The Hofman-Maldacena type limit: 



n 



L = 2, |J|>1, ^finite; A finite (2.15) 

where in the last two cases, we use a two-site lattice, corresponding to a two-magnon 
state. It is of course possible to study this in a more general L-site lattice. 
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2.2 Coherent String States 



We would like to be able to describe a coherent string state on the background given 
by an arbitrary distribution of black rings, as discussed in the previous section. In 
order to be able to do this we find a coherent string prescription which is normalizable 
on each of the m black rings separately. The coherent state satisfies the equation [23] 

a\z)i = z\z)i, (2.16) 

where the / labels the l^^ ring in which the wavefunction is well defined and normaliz- 
able. 

The coherent state can be expressed as 

k=l n=0 ^ k=l n=l ^ 2*^^ 

which is normalizible on the l^^ black ring, i.e. R2i~i < \z\ < R2i- The state |n, a^) 
corresponds to the n^^ excitation on the outer edge of the k*^ ring while \~n,hk) 
corresponds to the — n*'* excitation on the inner edge of the k^^ ring. It can be shown 
that the sums over k can be consistently restricted to the ranges in f l2.17p due to the 
normalizability on the l*^ ring. A more detailed derivation of the above is included in 
Appendix |Al 

In fl2.17p . we sum over all the positive occupancies, or in other words, positive 
angular momentum movers along the outer circles, outside the V"^ black ring, and 
over all the negative occupancies, or, negative angular momentum movers along the 
inner circles, inside the l*^ black ring. Note that such a summation is important 
for normalizability within the l^^ ring. The coefficients ff^~^ and ff^ correspond to 
the relative contribution to the wavefunction from the excitations associated with the 
{2k — I)*'* and 2k*^ circles. Equivalently these circles correspond to the hk-, au corners 
of the Young tableau from the gauge theory point of view. 

f l2.17p is a Fock space state which sums over all occupation numbers, as well as all 
types of oscillators specified by the Young tableaux corners bk,ak. The real, positive 
coefficients ff^'^Jf^ can be shown to satisfy the following relation 

m m 
k=l k=l 

where vl, is a mixing matrix which allows for mixing between different rings. The f's 
can then be written as 



2k-l\2 



^2k-l 



n 



l(-^2fc- 



/?2 



n 



p=li^2k- 



^Ip) 



(2.19) 

(2.20) 
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_ Q^n^;+i(<^afc ~ Cbp) YlpJi^f^jCg^ - CgJ (2 22) 

where we have used the identity (12. ip to write these expressions in terms of Cb^, Ca^,- It 
can be shown that these expressions are positive definite. By looking at the behavior 
of the norm {z\z) when merging two black rings, we get 

I m I m 

CRURlk-iUmk cnUCb.UCa, 

k — 2 k — / k — 2 k — / /o oo\ 

= —I - = — m - ' \^-^^) 
Y\Rlk n Rlk-1 NWCa.Yl a, 

k=l k=l+l k=l k=l+l 

where cr is an overall normalization which is unfixed from the coherent state condition 
but will be fixed later by putting a constraint on the integral over the inner product 
of such representations in the z-plane. Finally, one can write the coherent state wave- 
function on the entire z-plane as a sum over the states defined on each ring: 



m 



k) = El^)^ (2-24) 
1=1 

where 6x,y has unit support between x and y, and is otherwise. 

Having defined such a coherent state on the 1/2 BPS background, one can take the 
continuum limit: 

L>1, — <1; A<1, (2.25) 

as the product of coherent states on each lattice site: 

Gj \zj) = Zj \zj) . (2.26) 

This limit would correspond to the large angular momentum (L) limit on the string 
side. A general coherent state is then labelled by a coherent state on each site of the 
lattice: 

12:1,2:2, ...,zl) . (2.27) 

The continuum limit can be expressed in terms of a sigma model action with Hamil- 
tonian 02.71) . which is derived by looking at the large L limit of the lattice action: 



i=i 

■ /|2 



dt J da \^-{zck- zd,)\ogi{z\z)) - —\z'\' j , (2.28) 
V{z,z) = d^log{{z\z)) . (2.29) 
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This is a sigma model with target space the 2;-plane of the 1/2 BPS bubbhng chiral 
primary geometry [5], [6], see as well [7]. This Landau-Lifshitz type action coincides 
with the action from the string theory side for the string sigma model that describes 
a closed string traveling with large angular momentum L on the S'^ in the 1/2 BPS 
bubbling geometry under the uniform gauge, in which the angular momenta on the 
are uniformly distributed along the string. This uniform gauge corresponds to the 
uniform distribution of the Y fields on the lattice sites. This action can be written in 
terms of the above potential as: 



S 



2^ 

2^ 





(2.30) 




(2.31) 


-^(p'^ + pV^)), 


(2.32) 




(2.33) 



V = h'z, 

We have neglected higher order terms (-^) \z'\^, (^p-) which would correspond 

to the contribution from the two-loop dilatation operators and can be safely neglected 
in this limit. 

The equations of motion coming from this action are 

^z" ='-{d,V + d,V)z, (2.34) 

which have static solution 

= 0, p = 0, 0' = 0, p' = ci , (2.35) 
z = {cia + C3)e"'\ (2.36) 

where Ci,C2,C3 are constants of motion, and the solution corresponds to a straight 
string stretching across a black ring (and back again). 

One can also bring f l2.30p into the canonical form of a 2d complex scalar field theory. 
We can transform the z^if, bringing the kinetic term into its canonical form 

h^-{zz - zz) = -(^^ - ^^), (2.37) 

z = pe'^, ^ = re^"^, (2.38) 
hp = r, if = hz . (2.39) 

Under this transformation, the complex scalar field theory is described by the action 



L , 

S = — / dtda 
a' 



ivv - w) - T? fim )nW\ + ^(1^1 



2.40) 
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where 

fi\^\') = {drpf + ^, g{W?) = {drpY - 1 . (2.41) 

The potential terms in this action are dependent on the background in which the 
string is moving. This means that the action is dependent directly on the droplet 
configuration in the z-plane. In particular the topology of the background which is 
related to the number of inner and outer rings is captured in the potential. One can 
study the S-matrix of this model, via similar methods to [TB] for AdS^ x S^. For 
example for the AdS^ x case the potential is given by: 

- = |„2 \ 121 ' (2-42) 



- rr--2^"^^ ^ = re\ (2.43) 

z ^ if maps the disk into a whole plane and -Ri is the radius of the disk. For more 
general concentric black-ring configurations we have 



V 

z ^ 

k=l 



1 1 

+ 



R2 -\zr\ /?2 



(2.44) 



The limits of such an expression can be found in Appendix |Al 

We would now like to derive such a potential from the coherent state. We first 
consider the norm of the coherent states {z\z) . The coherent state fl2.17p defined on 
the ring is: 

I m 
k=l k=l 

where 



n,afc), (2.46) 



which have the property 



R2 _ U|2 p2 _ U|2 ' 

-^2k-l -"'2fc-l 



z,h\z,h) = - ,2 + 1 = - ^ ' , ,2 ^ (2-47) 



R2 |^|2 
-"-2fc _ H 

-^2fc ~ 1^1 ^2k ~ 1^1 



{z,ak\z,ak) = ,2 -1= ^0'^' , ,2 - (2-48) 



From the coherent state expansion we get 



\Z\Z 



\ _ ^Jl ' -'^2fc-l , V^Ui ) -"-2^ /'o /ln^ 



p2 _ ur ^ ' R2 

fc=l ^2k-l I-*! k=l^2k 
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which leads to a relation between the coefficients in the coherent state expansion 



2k\2 



(2.50) 



k=l 



k=l 



Note that the coherent state norm is nonzero only on the black rings. The expression 
f l2.49p can be summed and factorized in the numerator, and rewritten as 



i-i 



ciY[{Rik-\A') n (^2Vi - kr) (- kr) 

k=l k=l+l 



I 

r 

k=l 



(2.51) 



n {Rlk-i-\A')I[{Rlk-\A') 

k=l 



k=l 



n (1 



w — 



I 

n 

k=2 



r2 

-"-2fc-l 



m 

n 

k=l 



2\ ' 



(2.52) 



The norm of the coherent state can also be written in terms of the potential, V, as 



{z\z)i = e 

which can be used to calculate the potential: 

i-i 



fV{z,z)dz 



(2.53) 



V{z,z) 



k=l 



Z 



d2 

^2k-l 



R2 
^2k 



m 

E 

k=l+l 



P2 

^2k-l 



Rlk 



d2 
^21-1 



Rli 



(2.54) 



Note that the existence of cr is consistent with the existence of an integration constant 
in f l233|) . 

The coherent state norm is non-zero only on the black-rings as the potential diverges 
at their edges. This calculation of the potential matches with the V function derived 
in the sigma model from the string theory side. This shows that the coherent state 
wavefunction knows about the geometry of the multi-ring background. Moreover, the 
coherent state wavefunction has a different expression on each disconnected black ring, 
meaning that it knows about the disconnectedness of the droplets, i.e. topology, of the 
background. 

Another property we can study using the coherent state we have described above is 
the average value of the occupation number on each lattice site. From ( I2.17p . we have 
the coefficients in the coherent state expansion expressed as 



2k-l 



{z,bk\z)i 
{z,bk\z,bk) 



5 J I 



fl 



2k 



{z,ak\z)i 

{z,ak\z,ak)^ \z 



p2n 
^2k-l 



\2n 



2n 



r?2n 
^2k 



{z,ak\n,akf . 
(2.55) 
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The states are orthonormal, such that 

(— ni, bkl\ — n2, hk2) = 1 ■ ^ni,n2^6fci,fefe2' ('^l' (^klWli O-kl) = 1 ' '^ni ,n,2 (^afci ,0^2 • (2.56) 

Excitations of a different level or on a different ring are orthogonal. 

The coherent state is a superposition of states with different occupation numbers so 
in order to find the average occupancy we have to look at the probability of occupancy 
—n or n on a given inner or outer ring. For the inner rings the probability P{—n, bk) 
of occupation number —n on the fc*'* ring is given by: 



{z,bk\z,bk) 



P{-nM) = ' ] ,. = t¥ 1 - ' (2-57) 



and the average occupation number nfe^, on the inner edge of the k^'^ ring is (for n ^ 0) 

= l](-r^)P(-n, fe,) = • (2.58) 

n=0 1 

Similarly for the outer edge of the k^^ ring we have a probability P{n, ak) of occu- 
pation number n: 

I I \2 I |2(n-l) / I |2\ 

and the average occupation number on the outer edge of the k^^ ring is (for n > 0) 

OO ^ 

We can now define an occupation number operator n. The average occupation over 
all inner and outer rings (corresponding to the average U{1) J charge on the string) 
for the entire coherent state is 

/ OO m OO 

/ 1^1 \ E E (-'^>^fek)?(-^) + E E (^.«fck)?('^) 

/ \ _ \^ \ ^ \ ^li _ fc=in=o k=in=i ^na^\ 



E Uf'-y {z, bk\z, bu) nt, + Eifi'? ^k\z, au) n^, 

k=l k=l 



(2.62) 



\z\ 



where the norm can be written as 

I 



d,,. log {z\z), = \zf hi (2.63) 



{z\z), = Y^Ui'-'f h\z, bk) + Y^Uf^f "'^l^' ' (2-64) 



k=i 
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So we see that 

{n)i = \z\^ d\^\2\og{z\z)i , (2.65) 

which is again a quantity depending on both the position of all the black rings, and on 
the ring from which the measurement is being made. The average occupation number 
gives the average f/(l) J charge of the single trace part of the operator fl2.3p . 

{n) = ^, (2.66) 
and also gives the average angular momentum of the string on the 2;-plane. 



2.3 Density Matrix 

The states we consider in (12.31) are in general non-BPS states. They are dual to 
strings on the bubbling AdS backgrounds described by systems of concentric rings in 
the LLM plane. In the situation where we have only the 1/2-BPS backgrounds, all 
distinct representations (described by the Schur polynomials) form an orthonormal 
basis. Having broken supersymmetry with the inclusion of the trace operator, such 
orthonormality is however no longer guaranteed. Let us consider two very similar 
background geometries which have a large number of identical coincident black rings 
and a small region where the ring structures differ. The coherent state wavefunction 
will be able to distinguish such small differences. The radii, defining a given ring 
configuration (12.171) . can be viewed as the continuous parameters in the coherent state 
string wavefunction, and we can study the entire class of background geometries from 
the point of the view of this wavefunction. 

Having obtained the coherent string state on a general background geometry given 
by a Young tableau representation R, we can write the string state in this background 
as 

\z,R), (2.67) 

where R labels the background geometry, and is seen in the set of parameters in the 
coherent state wavefunction (the R^s). There are various constraints we may place 
on the representation but for now we study representations with fixed total number 
of boxes A and N total rows. We would like to define a reduced density matrix by 
studying the inner product between different string states in different backgrounds 
(representations R) and then integrating out the string degrees of freedom. This can 
be achieved through the following definition of the reduced density matrix over coherent 
states, 

p{R, = ^dzdz (z, R\z, ^) = (-^' R^ 271 \z\d\z\ , (2.68) 

where because of the 6 functions in the definition of {z)^ the only contributions to 
the integral come from the regions in the LLM plane where the black regions in the 
representations R and R overlap. The wavefunctions are normalized (using the inte- 
gration constant cr in (I2.23P ) such that the reduced density matrix element of one 
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representation with itself is unity, 



{z,R\z,R)2'K\z\d\z\ = l. (2.69) 

Given a particular constraint on the representations (a constraint on the total num- 
ber of boxes, rows and columns) we can define an entropy by looking at all possible 
inner products between representations within a superselection sector. The entropy of 
such a system is given by 

S = -Tr{p\ogp), (2.70) 

where p = ^p, Tr[p) = D, Tr{p) = 1. In other words we take the trace over the 
space spanned by the Young tableaux basis (including the additional Young tableaux 
plus black inner disk configurations). 

Looking at a fixed total number of boxes and a fixed number of rows A^, if we define 
the number of Young tableaux with no coincident radii as Di, and the total number 
of Young Young tableaux as D = Di + D2, the entropy will be between \og{Di) and 
log(D). We can also put a constraint on the density matrix in terms of the size of the 
off-diagonal components. The diagonals of the density matrix are equal, due to the 
normalization described above. We define the density matrix as the addition of the 
diagonal part and an off-diagonal part /, 

P = ^{lD>cD + f), (2.71) 

S = \og{D)-^Tr{{I + f)\og{I + f)), (2.72) 

where D is the total number of Young tableaux. This density matrix measures the 
similarities of different Young tableaux, and due to the contribution from the off- 
diagonal components, the entropy gets reduced from log(-D). If the Young tableaux are 
all identical then the entropy reduces to zero, and the system becomes a pure state. 
If on the other hand, the off-diagonal contributions are small, which is the case for 
M ^ iV (a sparse distribution of rings with many large gaps), then the entropy is close 
to \og{D); in that case the second term in fl2.72p can be approximated as —-^Tr{p). 
We analyze the dimension D of a constrained set of partitions in Appendix |Bl In 
general we constrain the Young tableaux to have fixed and fixed M = 

In order to calculate the reduced density matrix we must understand the inner 
product between two non-equivalent representations. Consider the two coherent string 
states in the backgrounds with representation R and R with m and fh total black rings 
respectively. The coherent states are defined as: 

1=1 tn=Ofc=l ^ ^ n=lk=l V 2fc/ J 

(2.73) 
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mloo/ /A \ oom 



[=1 1^*1=0 fc=i \ " / "=lfc=/ 



■2i-l'-'^2f 



(2.74) 

The inner product between the two representations can be written in terms of a contri- 
bution from the inner edges b and outer edges a separately as there are no cross-terms 
between a's and fe's: 

(^z, R\z, = (z, R\z, Rj + (z, R\z, Rj ^ . (2.75) 
Looking at the a type contributions we find the following set of summations: 

m rh m rh oo ^ \ ''^ f \ " 

z,R\z,r) = EEEEE/f/^'M^ \^] {n,a,\n,a,)en.,_,,^J^ 

m rh m rh oo / i i2 

\z 



1=1 i=l k=l n=l \"^2fc 



R^^ _ j ^2k"--2k 

and similarly for the h type contributions 

z, R\z, R 



b 

m m m m oo / r> \ 



^■21-1^^2 



1=1 k=l k=in=0 

m rh m rh oo / i i2 

2 



1=1 f=i k=l k=in=l \ ^2fc-l 1^1 



The calculation of such an inner product is complicated both by the large numbers 
of summations in the above terms, and also by the products implicit in the f's. Having 
calculated the inner products the integration must be implemented. The integration 
itself is not complicated but because there are very many terms we have performed all 
such computations using Mathematica. In order to find the full reduced density matrix 
it is also necessary to find all the possible Young tableaux for a given set of constraints. 
We also do this using Mathematica. The algorithmic procedure is outlined in Appendix 
Owhile an analytical estimate of the dimension of the reduced density matrix is given in 
Appendix [HI The dimension of the reduced density matrix for A = 40, M = = 9 is 
greater than 800 meaning that there are about a million non-equal entries in the density 
matrix. Such a computation is possible but, for representations with M and A^ much 
greater than 10, the combinatorics very quickly become unmanag eable, D > C(e^). 

In the following we plot some of the results for various M and A^, noting that for 
M and A^ of order 10, we can consider that we are approaching the large A^ regime. 
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For M = 6, = 6, A = 18, there are 29 distinct representations (the representa- 
tions are given in appendix [C]) . The reduced density matrix can be diagonahzed and 
the logarithm of the 29 eigenvalues calculated. The non-diagonalized reduced density 
matrix along with the logarithm of the eigenvalues are plotted in figure [21 Note that if 
there were no mixing, all eigenvalues would lie on the dashed line. Figures [3] and H] show 




Figure 2: The values of the 29 by 29 elements of the reduced density matrix for the M = 
= 6, A = 18 ensemble. On the right is the plot of the logarithm of the 29 eigenvalues 
of this matrix. The dashed line is the line that the eigenvalues would fall on if there were 
no mixing. The entropy of this particular ensemble is 2.50. The maximum entropy for an 
ensemble of 29 elements is log 29 = 3.37. 

the same plots for M = = 7, A = 24 for which there are 81 distinct representations. 

In figure \5\ we plot the entropy as a function of M and for A = . We see 
that for large M and A^ the entropy goes much slower than logD meaning that the 
off-diagonal elements are extremely important in describing the purity of the state. 
Indeed the constraint we are looking at M = A^ and A = is one in which we would 
expect a large number of overlaps in diagrams. This means that, in this case, the 
number of degrees of freedom is actually much less than the dimension of the density 
matrix. 

The limit of small A^, large M, however gives a result much closer to log D as can 
be seen in figure [6] where the points are the values of the entropy calculated for A^ = 3 
and the line is given by logD, showing an almost perfect match in this limit of no 
overlaps. Note that there is a subtlety in such a calculation which is dependent on 
how one regularizes the divergences in the reduced density matrix. See Appendix O 
for more details. 
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Figure 3: The logarithm of the 81 eigenvalues of the reduced density matrix for the ensemble 
of representations with M = N = 7, A = 24. The dashed line corresponds to the positions 
of the eigenvalues if there were zero mixing, i.e. — log 81. 



2.4 BMN String States 

Having studied the continuum limit, we turn to the two-site lattice in the BMN limit. 
Such states can be viewed as the two-magnon states of the string with large J charge. 
We can project onto states with fixed angular momentum J by performing the contour 
integral over the product of two coherent states as follows: 

= y dzz-^^^^^ \z) \z) , (2.78) 

however unlike in the one black-ring case [6], l^;) is a sum of contributions from each 
ring, so we have to sum over rings: 

m m 

\z)\z) = Y.Y.\')Az),dmi'), (2.79) 

1=1 l'=l 

where 6{l) is the theta function with unit support only on the l*^ ring (this is simply 
a rewriting of fl2.24p without referring to the ring radii). 

Since we are only looking at a single representation, R, (both \z) s are the same), 
the product of 9's means that the above simplifies to: 

m 

\z)\z) = Y,\z),\z),e{l). (2.80) 
1=1 
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Figure 4: The reduced density matrix for M = A'" = 7, A = 24. We see the line of 1/81 on 
the diagonal and a very large non-trivial mixing on the off-diagonals. 

We therefore start by looking at \z)i\z)i 

n / I oo /pN"- moo / \ n \ 

p=l s=0 \ ^ J ,=1 \^2p/ J 

Only terms which have a single negative power of z will contribute to the integrand. 
Therefore, for positive J this is only possible for terms of the form \—n) \n + J) , 
\n + J) \ —n), \n) \ J — n) and \J~n)\n). There will be no contributions from a product 
of two negative movers (two h types) only. 

Calculating the contribution from each / and summing we arrive at the final ex- 
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Figure 5: Plot of the entropy for the superselection sector A = This plot goes, for large 
M and N roughly as 5" ~ log(Af A^) ~ log(log(L')). This is a lot slower than the log(D) which 
would be the answer with minimal mixing between states. The conclusion therefore is that 
the non-trivial mixing is extremely important for M ~ iV and A = Such a conclusion 

will not be true for sparse representations. 

pression for the angular momentum eigenstate as: 

m ( / I oo m ^ 

1=1 y \k=l n.=0 p=l 

(m oo I \ 
E E E ff'f?-"R~.t'^Ri-. I- + I--' b,) 
k=l n=0 p=l / 

^ m J— I m \ "\ 

E E E fi'fi'R2kR2p' \n, a,) I J - n, a,) . (2.82) 

,k=l n=l p=l / J 

The analysis for the total negative J states is similar, except that we now pick up 
products of two b type states, and do not pick up the products of two a type states. 
We can now use fl2.82p as a general ansatz to write down an energy eigenstate. 
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Figure 6: Point plot for the entropy as a function of M with = 3 and A = along with 
the line plot for \og{D — 1) where D is the number of elements in the ensemble for a given 
M, this value also matches log(M/4) to extremely high precision. The conclusion from this 
is that for M ^ N there is little mixing and the entropy is maximal. 

Simplifying this expression we arrive at, 

m m ( / oo 
fc=l p=l L \n=0 

+ Fl':;r'R2t'^Rl-i \n + J, -n),,,) + Fn'^'R-^kR^;' |n, J - n),, j | , 

(2.83) 

where we write the two occupation numbers together into one bracket and leave 
the information about which ring is being excited as a subsript after the ket. The 
F's are calculated after applying an energy eigenstate constraint on the kets. A full 
explanation of this procedure can be found in Appendix [Dl In order to find the correct 
constraints on the F's we have to simplify the problem by looking at a decoupling 
limit which reduces us to the two site case already discussed extensively in [B] . The full 
problem of finding the F's in complete generality is outside the scope of the current 
paper. 

The decoupling limit which we choose is the one in which a single black ring, the 
l^^ is decoupled from the rest. In f l2.83p . we pick k = p = I, and since R21-1 = \/CbjN 
and R21 = \/ Ca, /N we can define 7 = CbjNi, and the state and algebra are very 
similar to the one black ring case [6] . From the general expression (12.831) , we get in the 
decoupling limit that 

° / \-f -^"^ °° / \^ 

n=— 00 ^ ^ n=l n=J ^ ^ 
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where Vn = \ J — n,n). 

The energy eigenstate equation is 

H 10) = 2X{a\ai + a\a2 - 4ai - ajas) |0) = E |0) (2.85) 

i.e. 

(2\{a\ax + a\a2 - a\a^ - ajaa) - ^) |0) = . (2.86) 
Applying this operation to our state we get: 

{0 \ — — 

Vn-\ + t'n+l) 
n=— oo ^ 

J-1 

+ E/„((2 + 27t;„-(l+7)(t;„_i + ^;„+i)) 

n=l 

+ E(r^) /„ ((1 + 27)t;„ - v/7(l + 7)K-i + ^^n+i)) 

' ^n=-oo ^ n=l n=J ^ 




(2.87) 



This equation generates recursion relations between triples of /„'s and we can solve 
the whole system iteratively starting from just two e.g. /o, /i- These recursive 
equations will be almost the same as [6], except that in fl2.87p we have E = ^^E 

instead of E. We thus have three unknowns, fQ,fi,E. /o is set to 1 via an overall 
normalization freedom while /i and E are fixed by demanding convergence of the series 
in the two asymptotic regimes n — ±oo. The calculation is discussed in detail in the 
appendix but the outline is the following: Three generating functions are derived which 
describe in closed form the recursion relations between /q and /j+i and /oo and, 
/_! and f-oo- Using these generating functionals we are able to write /_oo and /oo in 
as functions of E and /i. Convergence of the sums puts a constraint first on /i from 
the n — )■ oo regime and finally a constraint equation is found for E from the n — >■ — oo 
regime. More detailed discussions are in Appendix [Dl One finally gets that for large 
J states, the eigenvalues are 

Ca^ 2Tx^n^\ ( 2 + 47\ _ Ca, Sn^n^X 2 + 47 



_ 1 _ = ^ : 1 - , (2.88) 



^ The state ^2, ?T-i) is equivalent to the state under a cyclic permutation of the integers, 
e.g. to TT-i, [6], due to the cyclicity of the trace part of the operator (|2.3p . For a two-site 
system, one can write a general state as the sum of a symmetric state, i.e. one with the property 
\ni,n2) = |n2,ni), and an anti-symmetric state, i.e. one with the property |ni,n2) = — |n2,ni). The 
energy eigenstate equation (|2.86p is to be supplemented by the projection onto symmetric eigenstates. 
Thereby one should project out half of the eigenfunctions of the equation, associated with h = 2n + l, 
and only keep the other half with n = 2n [6] , and interpret the n as the mode number of the BMN 
state. Therefore finally (|2.88p . (|2.89p are in agreement with [B]. 



20 



where here n = 2n, and n is the oscillator mode number in the BMN state [15]. A 
similar procedure for large — J states produces similar recursive relations and yields 



These are strings excited along the exterior and interior edges of the black ring. The 
different coefficients in the energy spectrum reflect the different effective radii of cur- 
vature when we take the plane-wave limit along the geodesies near the different circles 
[1] . Thus, by measuring and comparing the effective units of energies of the states, one 
can know the exact location of these (short) string excitations in the bulk spacetime. 

3 Discussion 

In this paper we have studied the quantum string excitations on top of a class of axially 
symmetric bubbling geometries dual to the 1/2 BPS chiral primary operators. This 
work is a continuation of the research undertaken in [S] , [B] , [T] , [S] . On the gauge theory 
side, we have studied operators of the form of a product of a Schur polynomial Xr{^) 
and a single trace of Z's and impurities Y's. The Schur polynomial part corresponds to 
a general Young tableau which specifies a background geometry. The single trace part 
can be viewed as a bosonic lattice of Y's, and for each site there are, in general, 2m type 
of bosonic excitations, m of which have positive occupation numbers n, corresponding 
to the m inward corners of the Young tableau, at which point one can add n extra 
boxes there. The other m, with negative occupation numbers, —n, correspond to the 
m outward corners on the Young tableau, where one can subtract n boxes. 

Developing the work of [5],[6],|^ we have solved the problem of the coherent state 
for a general geometry with arbitrary concentric ring distributions (i.e. a general set of 
radii), and equivalently for a general Young tableau. We have also studied the algebra 
of the shift operators (or ladder operators) a, a\ acting on the different types of bosonic 
excitations on the lattice, and extended the algebra of [7] by including the mixing of the 
|0) states of different black rings. These different types of excitation correspond to the 
different possible circles in the geometry. The coherent state wavefunction has different 
expressions on the m different and disconnected black rings. On each black ring, the 
contribution to the Fock state wavefunction, normalizable on that ring, comes from all 
the positive mover states along the circles outside that ring, and all the negative mover 
states along the circles inside that ring. 

The norm of the coherent state wavefunction is related to the metric function V 
that appears in the dual string sigma model [5], [6], [7]. The radii that appear in the 
denominator of the norm are those radii for which the angular momentum excitations 
of the string are included in the coherent state expansion. The radii that appear in the 
numerator of the norm are all the rest of the radii in the geometry, and emerge due 
to the factorization of the numerator after summation. We have also calculated the 
average occupation number on each black ring, which are related to the metric function 




(2.89) 



21 



V. We have checked that the coherent state wavefunction behaves nicely under the 
merging of two nearby rings. 

Having obtained the coherent state wavefunction for a general concentric ring geom- 
etry, or equivalently, for a general Young tableau, we used the wavefunction as a means 
to study the different background geometries since the continuous parameters of these 
geometries are encoded in the coherent state wavefunction itself. We have then gone 
on to define a reduced density matrix element between any pair of geometries or Young 
tableaux. This has been implemented by integrating out the string degrees of freedom 
from the cross product between the coherent state wavefunctions living on these ge- 
ometries. The reduced density matrix measures the degree of similarity between any 
pair of geometries. In the most extreme case, if the two geometries are almost identical 
except for some very small inequivalent region, the matrix element between them is 
close to 1/D, where D is the total number of Young tableaux in a given constrained 
ensemble. We use computational methods to obtain both the representations in the 
ensemble and furthermore the reduced density matrix itself. This calculation is a novel 
development, which allows us to study the mixing of non-supersymmetric operators. 
The mixing of the tableaux is due to the non-BPS property of the entire operator 
(Schur polynomial times single trace). 

We have also studied an ansatz similar to the coherent states that describe the 
BMN states with large J charge. The coherent state has no fixed total J, but includes 
a superposition of all states with different J's. By picking the fixed J solutions from 
within the coherent state expansion, and using these states as the basis for the super- 
position to form a BMN eigenstate of the dilatation generator, we are able to solve 
the eigenvalue problem and find the eigenvalues. The spectrums encode the different 
radii of the rings within the geometry and therefore the locations of short string exci- 
tations in the bulk gravity dual. It therefore tells us some locality information of the 
bulk quantum gravity. We also show that one can obtain the BMN states with large 
negative J charge, which correspond to the negative movers along the inner circles of 
the different rings. For fixed large | J|, the number of sites of the lattice L corresponds 
to the number of magnons on the string. 

In finding the string sigma model in these geometric backgrounds we obtain a metric 
function, V , which describes the entire class of concentric ring geometries. The metric 
function ^ is a superposition of many terms each of which looks like a metric function 
for an AdS^ space with a different radius. A future goal would be to consider the 
solutions, and S-matrices of these sigma models in more detail. The straight string 
solutions are of particular interest, and may be related to the work of [21], Semi- 
classical long strings on these backgrounds were studied previously in [SB]- It would 
seem natural to consider these as operators of the form of the product of the Schur 
polynomial and an additional operator. 

A further direction of research would be to consider the higher loop dilatation 
generators acting on the lattice, having here extensively studied the application of 
the one loop dilatation generator. It would also be interesting to generalize the trace 
operator under consideration from one which contains only Z and Y impurities to one 
which contains fields from a larger sector, for instance the derivative operator. There 
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is still much to be done to fully understand many aspects of these types of operators. 
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A Derivation of the coherent states 

Here we give a more detailed explanation of the calculation of the coherent state, used 
extensively throughout this work. 

We assume that the coherent state on the Z*'* ring takes the form 



m 



k=l k=l n=l ^ ^ k=l n=l ^ ^'^^ 



where the excitations on the inner and outer rings are labelled by bk and and the 
ground state of the k^^ ring is labelled by |0)^. The ranges of the sum over k is due to 
the normalizability on R2i^i < < R2i- 

The algebras are written in section [2?T| and in particular, we assume that the action 
of the annihilation operators on a single positive mover gives a sum over zero modes: 

a|l,a,) = g^f^|0),, a|0), = ^|-l,&,) . (A.2) 

The mixing of the |0)^ states is natural, since when two nearby black rings merge, or 
equivalently, two nearby vertical edges of the Young tableau merge together, the |0) 
states on these two black rings, under this limit, have to be identified as the same state. 
The algebra without such mixing would reduce to m decoupled algebras, each of which 
has a similar form but with different coefficients. 

We then have that the action of the annihilation operator on the negative and 
positive movers respectively gives 

«f;(^) i-^'M=f:(^) ^\-n-iM) 

n=l ^ ' n=l ^ /V 

oo 



Z 

n=l 



E(^)V-,M-^|-1,M, (A.3) 
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n=l ^ ' n=2 ^ / V V 

oo ^ ^ n m 

= ^E^^ \n,a,) + zJ2vi\0)^. (A.4) 

, \ 2fc / ■_, 



The coherent state condition then needs to be satisfied which will provide a series of 
recursion relations for the f's. The condition 

«k)/-^k); = (A.5) 

yields the constraint 

m I mm m 

E^'^ 1-1, - E/""^ 1-1' + E/"(E-^- 10).) - E^^ 10). = • 

(A.6) 



fc=l ^ ^ fc=l fc=i J=l fe=l 



which gives a relationship between the f's and (7's 

9\ = iT~\ k^l; (7f = 0, A:>/. (A.7) 
We also have the following constraint coming from the sums over zero modes 

m m. I 

E/''(E^^|0).) = E/''"'|0), ' (A-8) 

fc=Z i=l i=l 

which can be rewritten as 

m m 

fc=/ fc=; 

and thus 

J - }!L - f2j-l_ \^ Jj_-j _ f2j-l -j _ f2j-l Jm~l f2j~l I'A in\ 

'^l ~ f2V '^l ~ Jl 2-^ f2k"k-- "m — Jm , ^^'m-l " f2m -I rn ' \^-^^) 

Taking the above constraints we are left with a simplified expression for the coherent 
state, given by: 

N), = E/SrT ( ^) i-".w + E/£E(i-) i".«'>. (a-") 

fc=l n=0 ^ ^ k=l n=l ^ 2^'^ 

where we have made a choice to absorb the zero modes in the hk states. From the 
coherent state expansion we get the norm of the state on the V"^ ring as 

I „2fc-lp2 m 2k p2 

k=l^2k-l 1^1 k=l^2k Pl 

= -(/^^-^)^ <t. = (/^^r, (A.13) 
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which can be re-summed and factored in terms of a quotient of products as 

l—l m 

n(^L-ki') n {Ri,.,-\zf)i-\zf) 

= c'-^^^ '-^^^^ ' (A-14) 

k=l k=l 

where q is an overall normalization constant. 

Since when there is no black region in the center, the norm is zero when z = 0: 

I m 
k=l k=l 

which has solution 

m m 

= -ci n (^2.-1 - Rip-X'''''"' n (^2.-1 - RipT''"'-'-'' (A.17) 

P=ly^k P=l 



rr2k\2 ^ (^lY[p=l+liRlk -^2p-l) np=l^fc(-^2fc Rip) ( tK ^R\ 

llp=ll^2fc - ^2p-lJ llp=Z^fcl^2fc ^2pJ 

m m 

= c,\{{Ri,-Ri_,)r'''-'' n (^2. - (a.i9) 

p=l p=l^k 

where sign(p — /) = 1, for p ^ I; and sign(p — /) = —1, for p < I. 
It can be seen that 

sign(cf_t"') = = (-l)#«>«2fe-i = (A.20) 

sign(cf^t^) = = (-l)#^>^2fe = 1. (A.21) 

The /f'^-^ and coefficients give the relative strengths of the response of the l^^ ring 
to the excitations on R2k-i,R2k circles in the m-black-ring geometry. 

It can be checked that they correctly produce all the non-zero coefficients of the 
2-circle case with 

ci 1 = 1 = Rl-Rf^ (A.22) 

since R^ — R\ is already proportional to N and we use the coordinate z in the units 
that Rl- Rl'ifil. 

It can also be checked that they correctly produce all the non-zero coefficients of 
the 3-circle case with 

p2 d2 p2 p2 

ci,2 = ^, c,, = ^. (A.23) 
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These match with the 3-circle result in [6|1 . Note our 4,3,2 is their 3,2,1, we have taken 
the Rl = hmit, and these normahzations make the norm of coherent state when 
Rf = 0, to be 1 (in the unit of cr) at the origin, z = 0. In this procedure, we first take 
Rl = 0, then z = 0. 

From this one can calculate the coefficient q, giving 

I m 

n^ifc n Rik-i 

k=l k=l+l 

where the subscript on cr denotes the representation R and not a radius. Before taking 
the Rl = limit the expression has the property that the coherent state wavefunction 
is normalized to zero at the origin, because there is no black region there. The limit 
from 2m circles to 2m — 1 circles is similar. 

The overall normalization factor for each black ring has the property that: 

(—1 m 

{z\z), = ^^f±L_ (A.25) 

U{Rlk-i-\A')U{Rlk-\A') 



(A.26) 

When one crosses from the ring to the (/ + 1)*'* ring, there is only a change of factor 

1 _ Al^\ (^ _ \A' 



Similarly, when one crosses from the l^^ ring to the (/ — 1)*'* ring, there is only a change 
of factor 

V ^ V _ ^^ 28) 

^21^1/ V ^21-2/ 

These are due to the sign changes of the contributions in \z) from the two circles that 
one passes by. 

The disappearance of the central circle leads to the limit 

1 , (A.29) 



(- 


\z 


?) 


{Rl- 




z\ 





'There is a slight typo in equation (58) of [6], where the prefactors of the last two terms are 
switched. 
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which gives a good hmit to the case of an odd number of circles - i.e. a disk plus annuli 
configuration. This is also a good limit for \z), since ~j^2Z^ ^ 0- 

A merging of two black rings when i?2j — > -R^j+i ^^^^ ^ natural and well behaved 
limit, as, e.g. 



|2 



^1, y (A-30) 



SO the formula reduces to the situation with one fewer black ring. This is also a good 
limit for \z), since c^^^^ — — c^"*, and contributions from the two merging circles cancel. 

The situation where the P'^ black ring is far away from others leads to the following 
limit 



(- 


\z\ 


?) 


(^1-1 - 


\z\' 




\z\ 





^21-1 ^21 1^ 



which reduces to the one black ring case, which is a form of decoupling limit for the 
jth ^--^^g 

The limit where we go from 2m circles to 2m — 1 circles takes the following form 

Rl^Q, V^V. (A.32) 

In fact the 2-circle case can be viewed as a limit of the 3-circle case with i?i = leading 
to the 1/z term. 

There is a further limit which would give the wavefunction in complimentary regions 
nomarlizable on the white regions instead of the black. Starting from the case of an 
even number (e.g. 2m) of circles, and sending the last circle to infinity, followed by 
switching the black and white regions leads to 

Rl^^oo, V^-V. (A.33) 

This procedure gives us coherent state wavefunction in white regions instead of black, 
as it is complimentary to the wavefunction in the black regions. 



B Dimension of the density matrix 

The dimension of the density matrix D in section 12.31 is related to the entropy of the 
superstar, S via 

= D. (A.34) 

The superstar geometry corresponds to a grey disk on the 2;-plane with a larger radius 
and lower fermion filling density than the AdS disk. It has the same N as the AdS 
geometry but has a larger radius on the 2;-plane |1] , [H] , [H] : 

2A 

Rlar = R\A^ + j^)=N + M, R\,s = N, (A.35) 

2A , 
M = — . (A.36) 
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A single superstar may be viewed as a microcanonical ensemble with total rows, and 
M total columns, both fixed. In other words, we have 

m m 

N = Y,N„ M = Y,M,. (A.37) 

i=i i=i 

The slope of such a Young tableau is, in the limit of large M and A^, 

M 2A , 

N 



and the filling density on the grey disk with the larger radius, is — jv+jw 

One can consider the entropy as the sum of entropies of each area element in the 
phase space, each of which is subject to a Bernoulli distribution {u, 1 — u}. The sum 
is written in [19] (see also |20]) which is 

S = --^ I d'^x\ulogu + {l-u)\og{l -u)] (A.39) 
ZTTn J 



= -N J d\z\^[u\ogu+ {l-u)\og{l-u)] , (A.40) 
where the area is quantized due to the quantization of the phase space 

d^x = M + N = N [d\z\\ = (A.41) 



2Txhj 7 ' ' ' 2^ 1/A^ 



see the footnote on page 3). For a uniform distribution in the phase space, m is a 

N 
M+N ' 



constant which is the filling fraction . ,^ ^ , and we can thus write the entropy as 



S = (M + N) 



M , fM + N\ N , [M + N 
log — ^7 — + „^ , ^. log 



M + N °\ M J M + N ° V A^ 
N + M 



(A.42) 



(A.43) 



where 



f ( ^ ] = + 
^ N + M J \ \ N Vm 



M , fM + N\ N , fM + N 
-log 77 — + ITT-^log 



M+N °\ M J M+N °\ N 

(A.44) 

This agrees exactly with [18], which is derived from the partition function approach. 
In the case of unrestricted A^ and M, the Hardy- Ramanuj an formula tells us that 

mm ^ 

W_v/A= WyMA^ = 2.5651 VA, (A.45) 
on the other hand from ( 1A.43I) 

g = /( /^ )v^^ 1.96052yA , (A.46) 
\N + M J ^ ' 
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which is consistent with the Hardy-Ramanujan bound. 
In the M ^ N regime (]A.43P has the hmit, 



Sr^N (^log^ + 1^ . (A.47) 

This formula is similar to the case of the type IIA geometry with units of D2-brane 
flux and K units of NS5-brane flux on the torus [22] . 

Another method to obtain this order of magnitude estimate for the entropy in the 
superstar geometry is by using Schur's theorem to find the number of partitions of A = 
with maximum rows. If D denotes the number of possible sets ni, rijv which 
satisfy the relation A = YliLi ^i^i non- negative, and {rj} are relatively prime, 

then 

[N — 1)! ri . . . 

in the large A limit. This should give an order of magnitude estimate, with the 
distinction coming from the fact that the set {rj} may not be relatively prime. If 
we let rk = k, which is the case we are interested in, we find 

logD- Ar(^logA + 2^ , (A.49) 

where we have used Stirling's approximation for the logarithm of the factorial in the 
large A^ limit. We see that ( IA.49P agrees with flA.47p in the leading order. 



C Procedure for computing the density matrix 

In this section we will give a detailed explanation as to how to calculate the reduced 
density matrix given a set of Young tableaux. We start with two representations shown 
in figure [3 The two ring configurations for the representations R and R are drawn 
only in their half planes so that we may see the coincident outer rings (marked with 
black lines) and coincident inner rings (marked with dashed blue lines). We may draw 
a more detailed picture of these two representations by simply plotting the radii in a 
linear fashion, as shown in figure [HI In this figure we have also labelled the and 
oscillators for the R representation along with the eight radii Ri,...,Rs and similarly 
the cifc, bk and Ri for the R representation. 

We see that the ring radii are coincident on two outer rings and one inner ring (the 
innermost). Note that the actual values of the radii are unimportant for this example. 
We simply aim to show the calculational procedure. 

We first write the coherent state wave function for the first ring structure (which 
has four rings) as: 

1^) = \z)^eR^^R^ + \z)20r.,^r^ + 12)3 6'^,,,% + \z)i9R^^R^ , (A. 50) 
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Figure 7: Two representations R and R given by their half plane disk diagrams. The black 
solid and blued dashed circles emphasize the coincident outer and inner rings. 



and for the second 



xy 



where each coherent state wavefunction \z)i has support only on the / ring 9 
for X < z < y and zero otherwise. Note that these are coherent states, though we can 
multiply the whole state by an arbitrary constant and it will still be a coherent state. 
We set this arbitrary normalization such that the inner product of the state with itself, 
integrated over the whole space gives where D is the total number of representations 
in the ensemble. 

J (fz{z\z) = j £z{z\z) = 1. (A.52) 
We can write out the coherent state wavefunction in full, in terms of the f's defined in 

( Km 

i ( 00 I \n oo4 



) = E EE/"-' ( ^ ) + EE/?' ( ■ 

l=\ Ln=Ofc=l ^ ' n=\k=l 



(A.53) 



4|ooi /A \" oo4 



i«> = E^EErM — ) i-">^> + EE/r|^) i"."^-)^"/^,-..^.- 

(A.54) 
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Figure 8: The two representations R and R compared in a purely radial diagram, showing 
their coincident inner and outer edges along with their a, b and Ri labels. 



Note that the fact that there are four rings in both representations is a coincidence - 
in many cases one looks at the overlap between representations with different numbers 
of rings. 

Now let's look at the inner product of these two objects. There is one clear simpli- 
fication we can make. This is that there are never any cross products between the a's 
and the 6's, so we can write: 

{^\^) = {^\z), + {z\z),. (A.55) 
Looking first at the a contribution in (1A.55P we see that we have the following sums: 

4 4 oo oo 4 4 



" 1=1 n=lm=l k=l ^ "^^^ \^2k) 

44oooo44 , xriZ-X"* 

EEEEEE^f^'M;^) ^^n^k,^RjR.^-r.RJk,,^,,k,, 



1=1 l^in=lm.=lk=l ^ ^'^^ \ 2fc, 

4 4oo44 _ /|^|\2n 

1=1 n=l k=l 



4 4 4 4 / I i2 

\z\ 



EEEE'^f'^'^^ d2 _ |.|2 ) ^R^-^RjR2i-i,RjR^i^^,R^j- 
1=1 1=1 k=l k=i ^ ^''i- 



^2k 



(A.56) 



There are two positions where the outer rings overlap, for / = / = 1 and 1 = 1 = 3. 
Note that it does not have to be the case that they are on the same rings. From this 
observation we get the following two terms contributing: 
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Z,R\Z,R) = fi fi I — —n 1 0Ri,R2^Ri,R2 

The first term includes a contribution from only one ring, wliereas tlie second term, 
wliicli comes from coincident a's on the third ring gives contributions from all the 
overlapping regions of all inward rings from that point. 

We can write a similar expression for the 6's and we find that the process is very 
similar except that the contribution comes from all rings out from the coincident radii. 
This time there is clearly no contribution from the second rings as they do not overlap 
at all, so there are again just three contributing terms here 

Z,R\z,R^^ = \_ (^flfl^RiR2^R^R2 + fsfs^RsRe^R^Re + flfl^RrRa^R^R^) ■ 

(A.58) 

Note that it is just a coincidence that the radii have the same label in the products 
of 6''s - they may come from overlaps of very different ring numbers in more complex 
examples. Because of the particular overlaps we see in the above diagram, we can 
simplify the 9, to get: 



z,R\z,r'^^ = |^ ^2 _ |^|2 j fifi^RuR2 + _ |^|2 j (jtfl,^RT.,R6 + fifi^Ri,R2^ 

(A.59) 

R\z, r)^ = (^^J^ j {flflOR^R,, + /3V3 ^i?si?e + flflOnrRs) ' (^-60) 

We will then need to take the integral over the complex z plane. The factors of 27r 
are not important because they will be taken care of when we normalize the diagonal 
of the reduced density matrix to where D is the number of diagrams contributing. 
The reduced density matrix therefore between these two terms, up to normalization is 
given by: 

(^RIR^^ = j d\z\ ^ ^J_^ \z\^ fifi^Ri,R2+ (^ j^2_!_ |_^|2 j (jtft^R5,R6+ fifi^RuR^ 

+ (^J^) (/i /i ^i^ii?. + flflOR.R. + flflOn.n}) ■ (A.61 

In fact this integral is divergent, but there is a very natural cutoff given by the fact 
that the phase space is discrete. The integration of the area in the phase space is really 
the sum of area elements in the phase space, i.e. d{\zY) = 5{R^). Given that the R^ 
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are quantized in units of j^, we subtract (or add) -^^ from the upper (lower) bound 
given by the i?2fc-i {R2k), giving: 



R\R 



flfl 



R2 



log [Rl 



Rl 



Rl 



+flfl\~^^og{Rl 



R5 



+ /T/iM-f i°g(^6 



+flfl (flog {R\ - \z\') + ^] ^ + flfl [flog (i?? - \z\') + 



+/4Vnf iog(i??-N^) + 4' 



K7 



This whole expression needs to be correctly normalized such that {R\R) and (^R\R 

are j^. Note that for small and M one should really apply a summation rather than 
an integration at this point. This distinction will not be important in the large M, N 
limit. There are however subtleties when we fix, for instance large M and small N due 
to coincident thin rings at very large radii. For small A^, these dominate when one uses 
the summation. If one considers an ensemble in which the total M is allowed to vary 
slightly, such outer ring coincidences are rare and this issue becomes unimportant. 

Having given an abstract example, we turn to a concrete calculation. In this case 
we look at the case for M = N = 6 and A = = 18. In this case there are a total 
of 29 Young tableaux which make up the ensemble. Some examples of these are shown 
in figure M The five diagrams with solid lines extending from the base are those which 
have central black disks within the annuli. We can calculate the reduced density matrix 
from these diagrams using the previously described algorithm and find the graph shown 
in figure [2j We concentrate on 2 points in the 29 by 29 matrix which correspond to the 
mixing between different representations. The two pairs of Young tableaux and disk 
configurations in the 2-plane are shown in figures [10] and [11] where the values of these 
reduced density matrix elements are -^0.12 and -^0.0014 respectively, where D = 29 is 
the dimension of the reduced density matrix. The fact that these off-diagonal values 
are non-zero shows that the entropy will be reduced from the maximal value of log(D). 
The mixing of representations via the breaking of supersymmetry leads directly to 
these non-zero off- diagonals. 

From the above discussion we can summarize the computational procedure for cal- 
culating the entropy of any chosen ensemble. The steps in the procedure are as follows: 

1. Choose the values of M, and A. Some of these may or may not be constrained 
depending on the limits of interest. 

2. Given the constraints on M, and A calculate all of the Young tableaux 
(including inner disk contributions) for these parameters. In general the number of 
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Figure 9: 16 examples of the 29 Young tableaux making up the ensemble of M = A'^ = 6, 
A = 18. 

diagrams, D, goes roughly as though the constraints on N and M can lower this 
number considerably. 

3. Look first at the diagonal elements of the reduced density matrix. There is an 
overall normalization of every \z) which is fixed by requiring that each diagonal element 
of the reduced density matrix should be j^. These normalizations are stored and are 
later used to weight the off-diagonal elements appropriately. 

4. Study the off-diagonal elements by cycling through all pairs of representations. 

5. Given two representations we first look for all coincident inner and outer edges 
which will contribute to the total reduced density matrix element. Calculate individu- 
ally the contributions from the inner and outer edges, integrating over the appropriate 
overlapping regions, summing to give the total and dividing by the weighting factors 
calculated from the diagonal elements. 

6. Calculate S = —Tr{p log p) by diagonalizing p. 

D Derivation of the BIVEN spectrum 

We include here an overview of the calculation of the BMN state using a series of 
recursion relations obtained from the energy eigenstate condition. The equation we 
discuss here is fl2.87p . This equation was obtained in |B] but the details of the calculation 
were not given there. We think that it might be useful to the reader to include more 
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Figure 10: The I'** and 4*'^ diagrams from figure [9] compared. The (1,4) reduced density 
element is -^0.12 showing small overlap in the representations. 




I I -3-2-1 12 3 I I I 



Figure 11: The I'** and 6*'^ diagrams from figure [9] compared. The (1,6) reduced density 
element is -^0.0014 showing almost no overlap in the representations. 

detailed intermediate steps in this appendix. Our results parallel theirs when the 
parameter is replaced by E. 

The constraint that the state is an eigenfunction of the Hamiltonian relates triples 
of f^s. We are able to solve the system iteratively by starting from two f^s which 
we chose to be /o and /i. The procedure has to be done in parts because different 
ranges of n have distinct recursion relations which cannot be solved using a generating 
function formalism. We start by looking at the modes from n = 2 to J — 3 (to be sure 
that we are not taking in any of the behaviour for the J*'* mode which has a slightly 
different recursion relation due to the factors of 7 coming both from the prefactor and 
the action with H. We also define ^E = E in fITHTD . 

We use a method of generating functions which can be defined from a general 
recursion relation. For the case where we have relations between triples of coefficients 
in our recursion relation we have an equation of the form: 

afn+i = bfn + cfn-1 , (A.63) 

and define the generating function which we would like to solve for as 

00 

B = Y,fnx\ (A.64) 

n=l 
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we can multiply the above recursion relation (]A.63P by x" and make the sum: 

^ oo oo 

- ('fn^'" = bB + xJ2 c/na;" , (A.65) 

leading to the relation: 



X 

n=2 n=0 



Solving for B we have: 



-{aB-afix) = hB + x{cB + cfQ). (A.66) 

X 



^ xcfn + afi / » „s 

B = — . (A.67) 

-a — b — cx 



In the case of the energy eigenfunction equation, concretely we have that: 

2x(7 + 1)a(x /^/o-/i) 

B = ^ , (A.68) 

2(7 + l)A(a; - 1)2 + Ex 

which we can write in the form of a summation by re-summing the geometric series, 
such that: 



n=0 



where 

4A(7 + 1) + V^a/^-8A(7 + 1) - E 

^- = «FTT^ ■ f^-™' 

By solving the recursion relations up to, for instance n = (J — 2) we can make sure 
that this gives the correct closed form for the recursion relation. To calculate /„ one 
simply calculates the coefficient of the x"' term. Now we calculate up to /j+i in 
terms of /o and /i using the recursion relation calculated and the generating function 
to take us to fj-2, and find that: 

r ^ r;^'^ (7 [rl - rl^) ((7^ + r+ (7^r+ - l)) /p + 7+r+/i)) 

1' {rl - 1) 

r-'-' (VttT (7 jrl' - rj) fp + r+ (rf - l) (7+ + r+ (7+r+ - l)) /Q) 

y {rl - 1) 

(A.71) 

where r_^_ = ^ and 7_,_ = 7 + 1. Now we perform the same procedure as before with a 
new generating functional, but going from n = J+2 to 00. This gives us the asymptotic 
behaviour for /„ which we use to fix fi. The generating functional, the equivalent of 
equation flA.GSp . is given by 

X{2x^-f+ + 7(2 - Ax) - 2x) + Ex 
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from which we calculate that the n*'*> J*^ term is given by: 

(7/2+) -'^ {R+1+) {rR+^' {ifj+1 - R+i+fj) li + I'^'R^;" if J - R+fj+i) r+) 



where 



Rll^ - 7 

(A.73) 



2A (7 + 7+) - M fx" - (7 + 7+)) +E^-E 

= -^-^x • (A-74) 

4A7+ 

We want the eigenstate to be normalizable which puts a constraint on the n — ?■ 00 
modes and leads to the constraint that: 

7/j+i-i?+7+/j = 0. (A.75) 

This is only possible for a particular value of /i and /q. We have the freedom, by the 
linearity of the eigenvalue equation, to set /o = 1 and so we are left with a constraint 
on fi which reads: 

^ ((^2J _ ^^^2 ^ (1^2 _ ^2jj {-f^R+ + 1) r_ + {rl^ - ri) (7+ - ^TTT)) 

" r-^ (-r_ {r'J - 1) (7+/?+ + l) + {r'J^' - l) 7+ + (r^^ " r'-) (7+ " v^)) 

(A.76) 

We can check that this is correct by performing the full recursion order by order in 
Mathematica up to a high value of n and checking that the series really does converge 
for fi of this form. Indeed the system is extremely sensitive to tiny changes in fi and 
the above solution gives precise convergence. 

Now we must perform the same procedure but working to n — )■ — cxd which leads to 
a constraint on E itself. Here using the recurrence relation for negative n's we get the 
following generating functional: 



n=0 



X (-7V+ - 7P+ + 7V+) {-if-ix^ - f-ix^ + 7/-2a^) 
7(7 + 1) (p+ - x) (xp+7 - 7 + xp+) - pI^ 



(A.77) 



where: 

47A + 2A - JaX'^ -8E-fX-4EX + E^ - E 

^ 4(.A.A) ■ (^^^«' 

From this we can pull out the n^^ coefficient in the series expansion to find the closed 
form value for /_„. Stepping through the n = —2 and n = —1 terms from the n = 
and n = 1 we can write this in terms of /i which has been fixed above. At this stage 
we have a closed form for /_„ in terms of E, J and n. This means that we can look 
at the asymptotics n ^ 00 and ask for a constraint on E such that the series is also 
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I rhs(E) 




Figure 12: Numerical calculation of the eigenvalues of the two site lattice Hamiltonian. 



convergent in this limit. Numerical studies show that the eigenvalues, E, coming from 
this constraint give precisely the same values as those from equation (32) in [6]. 

The numerical method is performed as follows. By taking the equation for E, 
similar to equation (32) in [6], and plotting the right hand side as a function of E and 
finding the points that it crosses we are able to find the eigenvalues. In figure [12] we 
plot / = a/— 1 times the right hand side of the equation versus E (solid line). The dots 
mark the points: 



J2 V J 

Ca, 2TT^n^X 2 + 47 



As discussed in section |2]4l one keeps only n = 2n states as the physical BMN states, in 
agreement with [6]. Figure [12] is plotted for J = 30, A = 2, 7 = ^. Note that one only 
expects agreement for the first few points for finite J due to the number of solutions 
to the transcendental equation. 
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